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A complete set of fermion and Higgs superflelds is introduced with well- 
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defined SO(W) properties and U(l) x Z2 x Z2 family charges from which the 
Higgs and Yukawa superpotentials are constructed. The structures derived 
for the four Dirac fermion and right-handed Major ana neutrino mass matrices 
coincide with those previously obtained from an effective operator approach. 
Ten mass matrix input parameters accurately yield the twenty masses and 
mixings of the quarks and leptons with the bimaximal atmospheric and solar 
neutrino vacuum solutions favored in this simplest version. 

PACS numbers: 12.15Ff, 12.10.Dm, 12.60.Jv, 14.60.Pq 



In a series of recent papers [ijj - || the authors have shown how fermion mass matrices can 
be constructed in an 5*0(10) supersymmetric grand unified framework by use of a minimal 
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Higgs structure which solves the doublet-triplet splitting problem ||. The construction 
was carried out in an effective operator approach with phenomenological input, including 
the Georgi-Jarlskog relations ||. Here we show how one can introduce a set of matter 
and Higgs S'O(IO) superfields with U(l) x Z2 x Z2 family charges from which the derived 
Higgs and Yukawa superpotentials uniquely give the structure of the fermion mass matrices 
previously obtained. The quark and lepton mass and mixing data are reproduced remarkably 
well with the solar neutrino vacuum solution preferred, provided the up quark mass is not 
zero at the GUT scale - otherwise the small angle MSW solution [0] is obtained. The right- 
handed Majorana neutrino matrix arises from a Higgs field which couples pairs of superheavy 
conjugate neutrino singlets. 

We begin with a listing in Table I. of the Higgs and matter superfields in the proposed 
model along with their 

Higgs Fields Needed to Solve the 2-3 Problem: 

45 B _ L : A(0)+- 

16: C'(l-p)++ 

16: £(-§)++, C\-\-p)-+ 

10: T x (l)++, T 2 (-l)+- 

1: X(0)++, P(p)+-, Z 1 (p)++, Z 2 {p)++ 

Additional Higgs Fields for the Mass Matrices: 

10: T (l+p)+-, ^(l + 2p)+-, 

T (-3+p)-+ ^(-l-3p)-+ 
1: F(2)-+, Y"(2)++, S{2 - 2p)-, 5'(2 - 3?)", 
Vm(4 + 2 P )++ 

Table I. Higgs superfields in the proposed model, 
family charges. As demonstrated in ||, in order to do all the symmetry breaking, one 45 
adjoint Higgs with its VEV pointing in the B — L direction, a pair of 16 + 16 spinor Higgs, 



plus a pair of 10 vector Higgs and several Higgs singlets are required. In order to complete 
the construction of the Dirac mass matrices, four more vector Higgs and four additional 
singlets are needed. Finally, one Higgs singlet is introduced to generate the right-handed 
Majorana mass matrix. 

From the Higgs SO (10) and family assignments, it is then possible to write down explic- 
itly the full Higgs superpotential, where we have written it as the sum of five terms: 

Wffiggs = W A + W CA + W 2/3 + W Hd + W R 

W a = trA 4 /M + M A trA 2 

W CA = X(CCf/M 2 c + F(X) 

+C'(PA/M 1 + Z X )C + C(PA/M 2 + Z 2 )C (\) 

W 2/3 = T X AT 2 + Y'Tl 

W Hd = TiOUY'/M + T CC + T (T S + T S r ) 
W R = T T' V M 

The Higgs singlets are all assumed to develop VEV's at the GUT scale. Wa fixes (A) 
through the Fa = condition where one solution is (A) oc B — L, the Dimopoulos-Wilczek 
solution ||. Wca gives a GUT-scale VEV to C and C by the Fx = condition and also 
couples the adjoint A to the spinors C, C, C and C without destabilizing the Dimopoulos- 
Wilczek solution or giving Goldstone modes. W 2 / 3 gives the doublet-triplet splitting by the 
Dimopoulos-Wilczek mechanism. Wh d mixes the (1,2,-1/2) doublet in 7\ with those in 
C (by F(j = 0), and in To and Tq (by F^ o =0). To fill out the model, we specify the 
50(10) x U(l) x Z 2 x Z 2 quantum numbers of the various matter fields in Table II. We 
require three chiral spinor fields 16j, one for each light family, two vector-like pairs of 16 + 16 
spinors which can get superheavy, a pair of superheavy 10 fields in the vector representation, 
and three pairs of superheavy 1 — l c fermion singlets. 
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l 6l (-|-2p) + - 162(-|+p)++ 16 3 (-i) ++ 
16(-I-p)-+ 16'(-|)- + 
16(±)+- T6'(-| + 2p)+~ 

10i(-l-p)- + 10 2 (-l+p) ++ 

li(2 + 2p)+" l 2 (2-p) ++ 1 3 (2) ++ 

l$(-2-2p)+" l c 2 (-2)+" ic(_ 2 _p)++ 

Table II. Matter superfields in the proposed model. 
In terms of these fermion fields and the Higgs fields previously introduced, one can then 

spell out all the terms in the Yukawa superpotential which follow from their SO (10) and 

U(l) x Z 2 x Z 2 assignments: 

W Yukawa = 16 3 • 16 3 • Tx + 16 2 • 16 • Ti + 16' • 16' • T x 
+I63 • I61 • Tq + 16 2 • 16i • T + 16 3 • 16 -A 
+16! • T6 7 • Y' + 16 • 16 • P + 16' • 16"' • S 

+I63 • 10 2 • C + 16 2 • 10i • C + 10i • 10 2 • Y (2) 
+I63 • 1 3 • C + 16 2 • 1 2 • C + 16! • lx • C 
+1 3 -1 C 3 -Z+1 2 -1 C 2 -P + 1 1 -1 C 1 -X 

where the coupling parameters have been suppressed. To obtain the GUT scale structure 
for the fermion mass matrix elements, all but the three chiral spinor fields in the first line 
of Table II. will be integrated out. The right-handed Majorana matrix elements will all 
be generated through the Majorana couplings of the V M Higgs field with conjugate singlet 
fermions as given above. 

With R-parity conserved, d — 4 proton decay operators are forbidden. The d — 5 proton 
decay operators induced by colored-Higgsino exchange that are generally present in unified 
models are present here but are not dangerous. It can be shown that the family charge 
assignments prevent any new and dangerous proton decay operators from arising. 



The procedure for deriving the Dirac mass matrices U, D, L, and N is the following. 
For each type of fermion /, where / = ul, u c l , d^, d c L , £l, f£, v L and z/£, the superheavy 
mass matrix connecting the / to the SU(3) x SU(2) x [/(l)-conjugate representation / is 
first found from Eq. (2) by setting the weak-scale VEV's and the intermediate-scale VEV, 
Vm, to zero. This will give three zero mass eigenstates for each type of /, corresponding to 
the three light families. Then the terms in Eq. (2) involving (Ti), (C), (To), and (Tq) give 
rise to the 3x3 Dirac mass matrices coupling ul to u c L , etc. This procedure is spelled out 
explicitly in |§. 

Under the assumption that the zero-mass states have their large components in the chiral 
representations 16i, I62 and I63, and all the other components are small, the Dirac mass 
matrices obtained have precisely the structure previously found in our studies by means of 



an effective operator approach: 
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with U and iV in units of Mjj and D and L in units of Mjj. The matrix parameters are 
identified with the Yukawa couplings and Higgs couplings and VEV's as follows: 

Mu = (^3)5(10)) M D = (£3)5(10), 

eMu = |3(a g /p)(t 2 )5(io)|, eM D = |3(a ff /p)(t 2 )s(io)|> 
V M V = (y'/s") 2 (t') 5m , aM D = -(c/y)(c')- 5{m) , (4) 

5M D = tato/s, 
6'M D = (t' i /s')e-^, 

where the subscripts on t 2 , h, t' and d refer to the SU(5)[SO(10)] representation content 
of the VEV's. The following shorthand notation has been introduced 



*3 = 


Al6 3 16 3 Ti(7l), 


*2 = 


A162I6T1 (Ti), 
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Al6'16'T! (Tl) , 
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Ai6 3 io 2 c'(C / )) 
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Al6 2 10ic(C)? 
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The parameter rj is introduced to give a tiny non-zero mass to the up quark at the Aq 
scale. Its appearance in N will also play an important role in the determination of the 
type of solar neutrino solution. It should also appear in D and L but its effect is negligibly 
small there and of no consequence, so it is dropped. The only phase then appearing in the 
matrices is associated with 5', as other phases are unphysical and can be rotated away 
with the exception of that associated with e. It turns out, however, that the best fits to the 
data prefer a real e. Hence <fi which can be identified with the complexity of the VEV of the 
S' Higgs singlet is solely responsible for CP-violation in the quark sector. The structures 
of the matrix elements given in Eqs. (3), (4) and (5) can be understood in terms of simple 
Froggatt-Nielsen diagrams [|H]] given in ||. 



Note that the 33 elements of the Dirac mass matrices are scaled by the VEV's of the 
10, T\. But the F = conditions for the Higgs superpotential require that the pair of Higgs 
doublets which remain light down to the electroweak scale arise from 5(Ti), 5(Ti), 5(C) 
and, to a very small extent from To and Tq terms, which are ignored here. In particular, we 
can write in terms of a mixing angle 7 

H v = h(T x ), Hd = 5(7i) cos 7 + 5(C) sin 7 (6) 

whereas the orthogonal combination has become superheavy at the GUT scale. Thus the 
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ratio of the 33 mass matrix elements found from Eqs. (4) and (6) is given in terms of the 
VEV's, v u and Vd of H v and H D , respectively, by 



Mu/M D = t> u /(f d cos7) = tan/5/ cos 7 



(7) 



Hence we find that the large My /Mr, ratio required for the top to bottom quark masses can 
be achieved with a moderate tan/3 provided cos 7 is small. 

Turning to the right-handed Majorana mass matrix, we use the zero mass left-handed 
conjugate states that were found implicitly above for the Dirac matrix iV to form the basis 
for Mr. The right-handed Majorana neutrino matrix is then obtained from the last two 
terms in Eq. (2), and we find 

/ At 3 0\ 
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where 

Ar 

Ae 3 A R = Xio 1 c Vm (V m ){ci/x)(c 2 /p22) 
Note that the whole right-handed Majorana mass matrix has been generated in this simple 
model by the one Majorana VEV coupling superheavy conjugate fermion singlets. By means 
of the seesaw formula ||11|| , one can then compute the light neutrino mass matrix 

/ 
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We now address the predictions of the mass matrices. For this purpose it is convenient 
to replace the parameters 5 and 5' by 
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which are essentially the left-handed and right-handed Cabibbo angles. In terms of the 
dimensionless parameters e, a, tL, £r, e ld , rj, A, and Mjj/Mrj, we then find at the GUT 
scale 
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Note that the Georgi-Jarlskog relations 0, = and = 3m", emerge as required 



Ae 2 
~ 1™Q 



(12) 



by design. The quark and charged lepton data are best fit at the low scale (see below) 
by assigning the following values to the model parameters: e = 0.145, a = 1.78, ti = 
0.236, t R = 0.205, 6 = 34° (corresponding to 5 = 0.0086, 5' = 0.0079, = 54°), 77 = 
8 x 10~ 6 , and M V /M D ~ 113. 

As noted earlier, in order to obtain the simple mass matrices in Eq. (3), we had to assume 
that the zero-mass states have their large components in the chiral representations I61, I62, 
and I63. The conditions on the state normalization factors are all satisfied provided the 
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following ratios are much less than unity: 

(a/p) 2 , (y'/s") 2 , (c/y) 2 , 

(13) 

(ci/x) 2 , (c 2 /p 22 ) 2 , (c 3 /z) 2 < 1 

With the numerical choice of parameters given above and near equality of the various Higgs 
couplings, we find {a/p) 2 ~ 0.02 and (y'/s") 2 ~ 6 x 10~ 6 , so the first two conditions are 
easily satisfied. Requiring that (c/y) 2 <C 1 and with the expression for a obtained from Eqs. 
(4), we find 

tan7 ^ §gf»a 

(14) 

tan/3 ~ \J o 2 + l(cos7)m°/m° <C m®/m b 

in terms of the T\ — C mixing angle, 7, in Eq. (6). With c/y = 0.1, for example, tan7 ~ 18 
which implies tan/3 ~ 6, a very reasonable mid-range value allowed by experiment. The 
others can also be satisfied ||. 

In |§ we have evolved the results in Eqs. (12) down to the low scales with a value for 
tan/? = 5, Aq = 2 x 10 16 GeV, A SUSY = m t (m t ), a s (M z ) = 0.118, a{M z ) = 1/127.9, 
and sin 2 6!jy = 0.2315. With the quantities m t (m t ) = 165 GeV, m T = 1.777 GeV, m M = 
105.7 MeV, m e = 0.511 MeV, m u = 4.5 MeV, V us = 0.220, V cb = 0.0395, and 5 CP = 64° 
used to determine the input parameters, Mu ~ 113 GeV, Mp ~ 1 GeV, and a, e, 6 



and r] given earlier, the following values are obtained compared with experiment |T2 
parentheses: 
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m c {m c ) = 


1.23 GeV 


(1.27 ±0.1 GeV) 


m b (m b ) = 


4.25 GeV 


(4.26 ±0.11 GeV) 


m s (l GeV) = 


148 MeV 


(175 ± 50 MeV) 


m d (l GeV) = 


7.9 MeV 


(8.9 ± 2.6 MeV) 


\V u b/V c b\ = 


0.080 


(0.090 ±0.008) 



(15) 



where finite SUSY loop corrections for m b and m s have been scaled to give m b (m b ) ~ 4.25 
GeV for tan/3 = 5. 



The effective light neutrino mass matrix of Eq. (10) leads to bimaximal mixing with 
a large angle solution for atmospheric neutrino oscillations JD1 and the "just-so" vacuum 



solution |TJ| involving two pseudo-Dirac neutrinos, if we set = 2.4 x 10 GeV and 



A = 0.05. We then find 

m 3 = 54.3 meV, m% = 59.6 /xeV, mi = 56.5 /xeV 
U e2 = 0.733, U e3 = 0.047, = -0.818, 5' CP = -0.2° 
Am 2 23 = 3.0 x 10~ 3 eV 2 , sin 2 26 atm = 0.89 (16) 
Am 2 2 = 3.6 x 10~ 10 eV 2 , sin 2 26 solar = 0.99 

Am 2 3 = 3.0 x 10~ 3 eV 2 , sin 2 29 reac = 0.009 

The effective scale of the right-handed Majorana mass contribution occurs two orders of 
magnitude lower than the SUSY GUT scale of Ac = 1.2 x 10 16 GeV. The effective two- 
component reactor mixing angle given above should be observable at a future neutrino 



factory, whereas the present limit from the CHOOZ experiment [15] is approximately 0.1 for 
the above Am| 3 . In principle, the parameter A appearing in Mr can also be complex, but 
we find that in no case does the leptonic CP-violating phase, 5' CP exceed 10° in magnitude. 
Hence the model predicts leptonic CP-violation will be unobservable. 

The vacuum solar solution depends critically on the appearance of the parameter rj in 
the matrix N, corresponding to the non-zero rj entry in U which gives the up quark a mass 
at the GUT scale. Should we set rj = 0, only the small-angle MSW solution [0 would be 
obtained for the solar neutrino oscillations. The large angle MSW solution is disfavored by 
the larger hierarchy, i.e., very small A value, required in Mr. 

In summary, we have constructed an explicit SO (10) supersymmetric grand unified model 
for the Higgs and Yukawa superpotentials which reproduces the fermion mass matrices 
previously obtained in an effective operator approach. All the quark and lepton mass and 
mixing data are fit remarkably well with a tan (3 in the range of 5 - 10 with matrix parameters 
which are also quite reasonable. 



10 



The research of SMB was supported in part by the Department of Energy under contract 
No. DE-FG02-91ER-40626. Fermilab is operated by Universities Research Association Inc. 
under contract No. DE-AC02-76CH03000 with the Department of Energy. 



11 



REFERENCES 



[1] C.H. Albright and S.M. Barr, Phys. Rev. D 58, 013002 (1998). 

[2] C.H. Albright, K.S. Babu, and S.M. Barr, Phys. Rev. Lett. 81, 1167 (1998); C.H. 
Albright, K.S. Babu, and S.M. Barr, Nucl. Phys. B (Proc. Suppl.) 77, 308 (1999). 

[3] C.H. Albright and S.M. Barr, Phys. Lett. B 452, 287 (1999). 

[4] C.H. Albright and S.M. Barr, Phys. Lett. B 461, 218 (1999). 

[5] S.M. Barr and S. Raby, Phys. Rev. Lett. 79, 4748 (1997). 

[6] H. Georgi and C. Jarlskog, Phys. Lett. B86, 297 (1979); A. Kusenko and R. Shrock, 
Phys. Rev. D49, 4962 (1994). 

[7] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); S.P Mikheyev and A. Yu. Smirnov, Yad. 
Fiz. 42, 1441 (1985), [Sov. J. Nucl. Phys. 42, 913 (1985)]. 

[8] S. Dimopoulos and F. Wilczek, report No. NSF-ITP-82-07 (1981), in The unity of 
fundamental interactions, Proceedings of the 19th Course of the International School of 
Subnuclear Physics, Erice, Italy, 1981, ed. A. Zichichi (Plenum Press, New York, 1983); 
K.S. Babu and S.M. Barr, Phys. Rev. D 48, 5354 (1993). 



[9] C.H. Albright and S.M. Barr, |hep-ph/ 0003251 



[10] CD. Froggatt and H.B. Nielsen, Nucl. Phys. B 147, 277 (1979). 

[11] M. Gell-Mann, P. Ramond, and R. Slansky, Report No. CALT-68-709: Supergravity 
(North Holland, Amsterdam, 1979); T. Yanagida, in Proceedings of the Workshop on 
Unified Theory and Baryon Number of the Universe, KEK, Japan, 1979 (unpublished). 

[12] Particle Data Group, C. Caso et al, Eur. Phys. J. C 3, 1 (1998); M. Bargiotti et al, 
|hep-ph/ 000T291 . 

[13] Y. Fukuda et al., Phys. Rev. Lett. 81, 1562 (1998); Y. Suzuki, in Proceedings of the 



12 



WIN-99 Workshop, Cape Town, 25 - 30 January 1999, to be published. 

[14] H. Fritzsch and Z.Z. Xing, Phys. Lett. B 372, 265 (1996); E. Torrente-Lujan, Phys. 
Lett. B 389, 557 (1996); V. Barger, P. Pakvasa, T.J. Weiler, and K. Whisnant, Phys. 
Lett. B 437, 107 (1998); A. Baltz, A.S. Goldhaber and M. Goldhaber, Phys. Rev. Lett. 
58, 5730 (1998); H. Georgi and S.L. Glashow, |hep-ph/9808291 ; C. Giunti, Phys. Rev. 
D 59, 077301 (1999). 

[15] CHOOZ Collab. (M. Apollonio et al.), Phys. Lett. B 420, 397 (1998). 



13 



